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Abstract 



*~ 1 ' We calculate the leading-twist O(a%f3o) corrections to the B —} -k transition form factor /+(0) 

in light-cone sum rules. We find that, as expected, there is a cancellation between the O(a1j3o) 
corrections to fsf+(0) and the large corresponding corrections to Jb, calculated in QCD sum 
rules. This suggests the insensitivity of the form factors calculated in the light-cone sum rules 
approach to this source of radiative corrections. We further obtain an improved determination 

,y- . of the CKM matrix element \V u b\, using latest results from BaBar and Belle for /+(0)|V^6|. 
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1 Introduction 

In the last decade we have witnessed major advances in the efforts to overconstrain the sides of 
the unitarity triangle, in order to test the CKM (Cabibbo-Kobayashi-Maskawa) mechanism of 
the Standard Model (SM). However, one side of the common parameterisation of this triangle 
is given by |V^,|/|V^,|, where Vij are elements of the CKM matrix, and recent determinations of 
\V u b\ have uncertainties of approximately 10% [I], as opposed to the error on measurements of 
\V c b\ from the inclusive channel B — > X c lv which is below 2% 0. Since the inclusive channel 
b — > ulv is dominated by the large b — > civ background, a competitive determination of \V u i\, 
promising both theoretically and experimentally, is found via the exclusive semi-leptonic decay 
B — > irlv. This requires information about the relevant hadronic matrix element, parameterised 
by the form factors f + (q 2 ) and /_(g 2 ), 

(7r(p)\u^b\B(p B )) = (p B +p)vf+(q 2 ) + (PB-P)nf-(Q 2 ), (1) 

where pb and p are the momenta of the B and n mesons respectively and q 2 = (pb — p) 2 ■ The 
beauty of this channel lies in the fact that in the limit of massless leptons, applicable to I = e 
and p, only f+(q 2 ) is required [3], 
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where Gf is the Fermi coupling constant and \{q 2 ) = (rn 2 B + m 2 — q 2 ) 2 — Am^rn 2 . for masses 
rriB and m^ of the B and n mesons respectively. Therefore the extraction of \V u b\ relies on 
the theoretical prediction for a single hadronic quantity f + (q 2 ), possible via non-perturbative 
techniques such as Lattice quantum chromodynamics (QCD) (see e.g. Refs. [HIS]) or QCD sum 
rules on the light-cone (LCSR). 

Theoretical predictions are usually confined to a particular region of q 2 , for example LCSR 
are restricted to large recoil energies of the pion, corresponding to q 2 < 6 — 7GeV 2 , and Lattice 
results to small values of the pion momentuno, i.e. q 2 > 15 GeV . Experimentally the q 2 
distribution has been measured with increasing accuracy at CLEO 015], BaBar [9] [T0 | [TT 1, [T2 t [T| 
and Belle [131 1 14] . In order to maximally exploit these theoretical and experimental results, 
one requires a well motivated parameterisation for the q 2 dependence of /+(g 2 ). There are 
a number of approaches, either simple pole-type parameterisations as in Refs. |15t 116] . using 
dispersive bounds to constrain the coefficients of a series expansion as in Refs. |1T|. [T5] or using 
the Omnes representation as in Refs. [19^ I20|. In all these, the normalisation provided by the 
LCSR prediction at q 2 = OGeV 2 plays a crucial role. In fact, one can obtain \V u b\ directly from 
the model independent result for |V^fe|/+(0), first calculated in Ref. [21] by fitting such shape 
parameterisations to BaBar data [10] . 

Light-cone sum rules are an adaptation of the traditional QCD sum rules approach [22] [23], 
considering instead the correlator of the T product of two quark currents sandwiched between 
a final on-shell meson and the vacuum |24[ 125] . This can be expanded about the light-cone, 
in terms of perturbatively calculable hard scattering kernels convoluted with non-perturbative, 
universal light-cone distribution amplitudes. The correlator can also be expressed as the sum 
over excited states, the first being the B meson which is then followed by a continuum of states. 



1 Note that the form factor at q 2 — GeV 2 was recently obtained in a quenched calculation on a very fine 
lattice [6]. 



Then assuming quark hadron duality above a certain continuum threshold, one can subtract 
this continuum contribution from both sides. Borel transforming this relation then ensures that 
this assumption, and the truncation of the series, have a minimal effect on the resulting sum 
rule. 

We are interested in calculating the subset of two- loop radiative corrections to /+(0) pro- 
portional to 0o, assuming, as discussed in Sec. 13.11 that this is a good approximation to the 
complete next-to-next-to-leading order (NNLO) result. In addition to allowing an improved 
determination of \V u b\, our calculation will enable us to investigate the size of these radiative 
corrections in view of the sizeable two-loop contribution to f B in QCD sum rules J26H2I]. The 
magnitude of this contribution is thought to be due to coulombic corrections, as explored in 
e.g. Ref. [281 . The LCSR approach to form factors involves taking the ratio of f B f+(q 2 ), also 
affected by such coulombic corrections, to f B - We therefore test the argument that radiative 
corrections should cancel in this ratio, provided both quantities are calculated in sum rules. 

The current status of the LCSR calculation of f+{q 2 ) is as follows. The next-to-leading 
order (NLO) twist-2 corrections to f+(q 2 ) were first calculated in LCSR in Ref. [29] and the 
leading order (LO) corrections up to twist-4 were calculated in Ref. |30j . Since the LO twist-3 
contribution was found to be large, further improvements were made by calculating the smaller 
NLO corrections [16] . A more recent update where the MS mass is used in place of the pole 
mass for mj can be found in Ref. [3TJ [32] . 

The following paper is structured as follows: in Sec. [2] we introduce the necessary notation 
and establish the framework required for the calculation, including the expression for the one- 
loop correction at leading- twist; in Sec. [3] we present details of the two-loop calculation and 
describe the structure of the divergences of the bare result and the renormalisation procedure; 
a detailed analysis of our numerical results, with predictions for \V u i,\, can be found in Sec. HI 
finally we summarise in Sec. El 

2 Set-up of the calculation 

Such as to briefly introduce the LCSR approach to the calculation of f+(q 2 ), and the notation 
which will later be required, we consider the correlator of two quark currents sandwiched between 
the vacuum and pion, 

n M = im b f d D xe- l P BX {7r(p)\T{u(0h fl b(0)b(x)i l5 d(x)}\0), (3) 

= ( PB + p)^Ii+{p 2 B , q 2 ) + (PB ~ p)^-{p 2 B , q 2 ). (4) 

In the region around the pole at p 2 B = m 2 B , H + {j> 2 B ,q 2 ) can be expressed in terms of f+{q 2 ) and 
the B meson decay constant f B , where 

m b (0\di 75 b\B) = m 2 B f B . (5) 

Above the B meson pole the contribution of the hadronic states can be described by the spectral 
density phadi leading to an expression for the correlator of the form 

IMPW) = fBm 2 B J±^ + / ds^. (6) 

m B PB Js>m 2 „ s P B 



Alternatively, in the Euclidean region where p 2 B — m B is large and negative, using a light-cone 
expansion about x 2 = 0, the correlator can be collinearly factorised into perturbatively calcu- 
lable hard kernels Tj_ (u, /x 2 ) and non-perturbative light-cone distribution amplitudes (DAs) 
4>( n > (u, /jL 2 ) for a given twist n, via 

U + (p 2 B ,q 2 ) = J^ f duT + ^(u,p 2 B ,q 2 ^ 2 )^ n \u,p 2 ), (7) 

where u is the momentum fraction of the quark in the pion, and /i is the factorisation or 
renormalisation scale. This factorisation theorem is not proved to all orders, but can be verified 
at a given order in twist or perturbation theory by the cancellation of IR and soft divergences, 
the latter arising when the convolution does not converge at the endpoints. The leading-twist 
pion distribution amplitude, (j)(u,fj, 2 ), contains the distribution of the momentum fraction u in 
the pion's infinite momentum frame for the lowest Fock state. We postpone the discussion of 
DAs to Sec. 13. 3| and here simply state the definition, in the Fock-Schwinger or light-cone gauge, 
to be 

{ir(p) 1^2(0)7^75 d(x)|0) = -if^ [ due^^u, fi 2 ) + ..., (8) 

Jo 

where f n is the decay constant of the pion, u = 1 — u is the momentum fraction of the antiquark, 

and the ellipsis indicates the contributions at higher- twist. Making the substitution u = (m 2 — 

q 2 )/(s — q 2 ) in the leading twist contribution to Eq. ([7|), and taking the imaginary part, we can 

define the spectral density /?T2 at twist-2, 

f °° n 

n + (p 2 B ,q 2 )= / ds-^ Y + ..., (9) 

Jo s-Pb 

where again the ellipsis indicates the contributions at higher- twist. Equating the expressions for 
H+(p B ,q 2 ) in Eqs. © and © results in 

f 2 /+(<? ) f i Phad f°° r PT2 , 1fl , 

fB m B~2 12- + / „ ds ~ -JT = / ds — —JT- ( 10 ) 
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m B~PB Js>m 2 B s-P B Jo s ~Pb 



Above the continuum threshold so , a continuum of states contributes and the approximation of 
quark-hadron duality is thought to be reasonable, such that 

Phad = PT2<8>{s- s ). (11) 

Subtracting the continuum contribution and Borel transforming both sides results in the sum 
rule for f + (q 2 ), 

U(q 2 ) = j 1 ^ fi^e-^B)/^, (12) 

JBm B J m 2 

where M 2 is the Borel parameter. The uncertainty introduced in making the quark-hadron 
duality approximation is reduced by Borel transforming, and further by choosing sq and M 2 
appropriately such that the result for f+(q 2 ) is flat with respect to these parameters. 

Returning to the original definition of the correlator in Eq. ([3]), we consider the NLO cor- 
rections to the leading- twist term in the expansion about the light-cone x 2 = 0, calculated in 
Ref. |29| . In analogy to Eq. ([7]), we express the correlator in the collinearly factorised form, 

n>|,g 2 ) =J2 J duT^(u,v 2 )^ n \u,p 2 ). (13) 



We perturbatively expand the leading- twist contribution to the correlator, 

nj 2 = /d»72 2 W 2 )«K«>/» 2 ) (") 



n ? , + £nw + g) 2 iv / n«..., (i 5 ) 



where the tree-level term LT^ is 



-\Um b [ l du<t>(u,^)tr{ ^ _% +mb J }. (16) 

\JpB — UP) ~ m h 





Although the 0(a s ) radiative corrections to the correlator, involving six further diagrams, were 
calculated in Ref. [29], we include the following expressions here as they will be useful in pre- 
senting the NNLO results, 

Fj, (17) 
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where the normalisation W is defined as 
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C F U m b , 





(18) 

for Cf = 4/3. Fj contains the total contribution of the traces and fermionic propagators for the 
weak vertex correction, B vertex correction, box, b quark self-energy and light quark self-energy 
diagrams. We factorise the gluon propagator out of F^ so that our notation can be adapted to 
the NNLO calculation more easily. Defining F^ to be 

Fj = F, WV + < V + F BX + F SE + F LS E) (19) 

the contribution of individual diagrams in Feynman gauge can be expressed as 

F WV = tr[ t-u* i-k ! + up' + m b ^ B -ui> + m h 

M a (k — up) 2 M (q — k + up) 2 — m 2 (pb — up) 2 — m 2 

F BV = tr{l JB-uj + vfib -i>B-$ + uj>-m h k^-uj 

M M (pb — up) 2 — m 2 a (pb + k — up) 2 — m 2 (k — up) 2 

F BX = t r uj>-$ fa-up'-k' + rn, l + v$_ 

M a {up — k) 2 M (pb — up — k) 2 — m 2 (k + up) 2 



pS e _ f r j>B-uj> + m b -p'b + u^ + ^t-m b p'b ~ uj> + m b 
M ** (pb — up) 2 — m 2 a (pb — up — k) 2 — m 2 (ps — up) 2 — m 2 



As in previous calculations, we work in the limit that the light quarks are massless, i.e. p 2 = 0. 
Therefore -F^ SE , the contribution of the self-energy diagrams for the external light quarks, 
vanishes as discussed in Sec. 13.21 In this paper, to avoid repeating what already exists in 
the literature, we will only concentrate on the technical details for the O(a 2 fio) corrections. 
Details of the NLO and higher twist contributions incorporated into our numerical analysis are 
as given explicitly in Ref. |16j . 



Pb / 




Figure 1: Feynman diagrams for 0(a 2 /3o) corrections to IT^ 2 . From left to right, the B vertex 
correction, weak vertex correction, box and b quark self-energy diagrams are shown. The external 
quarks are on-shell with momenta as indicated and the dashed line represents the B meson. 



3 Radiative corrections at order ol^/3q 

3.1 Calculation of the fermion bubble diagrams 

In analogy to QED, where the running of the /^-function is connected to the photon polarisation, 
Brodsky, Lepage and Mackenzie had the idea of associating the running of the QCD /3-function 
with fermion loop insertions in the lowest order corrections [33]. The scale for a given process 
can then be set by demanding that this contribution to the two-loop corrections vanishes, a 
procedure known as BLM scale setting. Physically, such a renormalisation scale reflects the 
mean virtuality of the gluon propagator [34] . 

In Ref. [35], the technique of naive non-abelianisation (NNA) was proposed, where the com- 
plete NNLO result is approximated by calculating fermion loop insertions, as for BLM scale 
setting, and replacing Nj by its non-abelian counterpart — (3/2)/3o- This idea was supported 
by the observation that in a number of cases where the remaining part of the two-loop correc- 
tions could be calculated e.g. higher order corrections to observables from hadronic vacuum 
polarisation and to the pole mass, it was found to be small in comparison to the O(a s (3o) con- 
tribution 2 ]. Using the NNA technique, we calculate the O(a 2 /3o) twist- 2 contribution to /+(0), 
keeping in mind that the NLO corrections to the higher twist contributions have been found to 
be comparatively smal]j. The expression to be calculated takes the form, 



n ( 2 ) 



Mi dm 



(U,/i 



d D k r(e)r(2-e) 2 / -k 2 
(2tt) d r(4 - 2e) \4^ 



1 



a/3 



k a l 



k 2 



F 



(24) 



where FT is as defined in Eq. (|19[) . The relevant Feynman diagrams are shown in Fig. [TJ 

The calculation is similar to the one-loop case, however, the additional fermion loop induces 
two important changes. Firstly, the tensor structure of the gluon propagator changes from the 
form 

-i ( „„ k a k?~ 



-ig a ^ 



k 2 



a/3 



k 2 



(25) 



2 Further, in Refs. [361 I37j . this idea was used to extend the BLM scale setting, by resumming fermion loop 
insertions in the lowest order corrections to all orders. 

3 Note that the various contributions to / + (0) were studied in Ref. [31] in the pole and MS schemes, and while 
at LO the twist-3 are comparable to the LO twist-2 contributions (~ 40 — 50%), at NLO, in comparison to the 
twist- 2 (~ 10 — 20%), the twist-3 contributions are better under control (~ 2 — 4%). 



resulting in additional terms in the trace (although these cancel in the sum of all diagrams 
due to gauge invariance [38], serving as an additional check of the calculation). Secondly, the 
factor r(e) means that the integrals must be expanded to a higher order in e. The increased 
complexity of the calculation is slightly compensated by the fact that we set q 2 = 0, however two 
scales (p 2 B and m^) and one dimensionless parameter (u) remain. We perform the traces using 
the package FeynCalc[39], and expand the hypergeometric functions using the Mathematica 
package HypExp [3D]. The resulting analytic expression must then be simplified and rearranged 
into a form facilitating the convolution with the distribution amplitude. 

3.2 Structure of the divergences 

The bare 0(a 2 Nf) results for njt , contain both infra-red (IR) and ultra-violet (UV) divergences. 
These are treated in naive dimensional regularisation (NDR), with totally anti-commuting 75 
due to the presence of two 75 matrices in the trace, renormalising the UV divergences in the MS 
scheme. As mentioned earlier, in NDR the light quark self energy diagrams vanish, as the UV 
and IR divergences arising from these diagrams cancel. On adding all the diagrams together, 
we first perform the gluon self-energy renormalisation using the 0(a s Nf) contribution, ^3y M , 
to the corresponding renormalisation constant ZyfM |41j . 



4ym = ~Cf ( ^ ) , (26) 



2 

Je 

multiplied by n^ . The left-over UV poles are completely removed by mass renormalisation, 
using the 0(a 2 Nj) contribution, Z„ L , to the renormalisation constants Z m , 

Z ™ = Cf ("7 + I) ' (27) 

multiplied by 11^ . Collecting what we assume to be the remaining IR divergences in 11^' IR 
and subtracting this quantity, 

n ( 2 ),en. = n (2) _ z (l) Mn (l) _ fflrfO) _ rf2),T m (2g) 

leaves 11^ ' UV and IR finite, however we are still to determine the origin of the IR divergences 
contained in 11^ IR . 

3.3 Convolution and scale dependence 

The leading-twist pion DA defined in Eq. ([8]) can be expanded in a series of Gegenbauer poly- 
nomials, 

oo 

4>{u, v 2 ) = 6n(l - u) Y, a n (pi 2 )Cll 2 {2u - 1). (29) 

n=0 

Here a n are known as Gegenbauer moments, and in the case of the pion the odd moments are 
zero by G-parity. The expansion is usually truncated, as the higher moments are suppressed due 
to the highly oscillatory behaviour of the Gegenbauer polynomials. However, the truncation is 
only justified if the hard scattering kernel 7/1 is slowly varying and non-singular for all u [42] . 
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— -jCf + —Ca 


■i^F [V{n 1 -J 1 T-B 4 2(n+l)(n+2)J 



Table 1: One- loop anomalous dimensions of the parameters a n , 773,4 and ^3,4 describing the 

DAs [mi 



We include terms for n < 4 up to 0(a s ), but we assume that at 0(a 2 Nf) the effect of 02, 4(7/) is 
negligible^, and adopt the asymptotic DA (i.e. 0(u, 00) = 6u(l — «)) to simplify the convolution. 
As the previously calculated twist-3 and 4 contributions are included in our numerical analy- 
sis, the corresponding DAs are also required, as defined in Ref. [IB]. In the same reference it was 
shown that, for a given twist, the two and three particle distribution amplitudes can be related 
by an equation of motion, resulting in a reduced number of independent parameters: 773,4 and 
W3,4. These parameters, as well as the moments a n , are known to renormalise multiplicatively 
to leading log accuracy [42], 

,2\\ lc/Po 



c(M ) = c(/J ) 



a s (/i 



(30) 



where /Jo is the initial scale at which the parameter was calculated and 7 C are the one-loop 
anomalous dimensions defined in Tab. Q]for c = a n , 7/3,4 or 0Jz,A- 

Coming back to the renormalisation of our NNLO result, the UV structure of the asymptotic 
DA can be factorised into the function ZAu, v) [38]. This can be related to V(u, v), the evolution 
kernel governing the renormalisation group (RG) running of the asymptotic DA, via 



V(u,v) 



Z,p(u,v) V dfi 







(i- 



Z^,(u,v) 



(31) 



V(u, v) is defined in Refs. 
given to 0{a 2 s Nf) by, 



, where it was first calculated to two-loop accuracy, and is 



a.\ 2 l 



V(u,v) = ^-V (u,v) + ^) - Nf C F V N (u,v) + 



(32) 



Explicit expressions for Vo(u, v) and Vjv(u, v) can be found in Ref. [H], and the ellipsis indicates 
other 0(a 2 ) and higher order terms. Z, (u,v), i.e. the 0(a 2 Nf) contribution to Z^u^v), can 
then be reconstructed from the evolution kernel, and expressed in terms of Vq(u, v) and Vn(u, v), 



a 1 
Z^u, v) = 5(u, v) + - 1 - 2V (u, v) + 
47r e 



(^£\ 2 — 

UJ e 2 



N f C F (\v (u, v) + e V N (u, v) J + . . . . (33) 



(2,0), 



On convolution with the tree-level hard scattering kernel 7]} ' (u, /J ), i.e. the leading contribu- 
tion to Tji (u,fi 2 ) in Eq.(TT5]), the divergence up to 0{a 2 Nf) takes the form 

jtfUuv = f du f dv-C F V N {u,v)T^ Q) {u,ii 2 )(j){v,ii 2 ). (34) 



4 This can be inferred from Fig. 1 of Ref. [3], where the respective size of different contributions in a n to f+{q 2 ) 
were shown as a function of q . 



Note that the terms in Vq(u,v) are symmetric in u, v, and therefore vanish since we use the 
asymptotic DA. The UV divergence of the DA cancels the IR divergence of the hard scattering 
kernel exactly at 0(a 2 Nf), i.e. Yi^ = — Txjj, IR . Therefore the IR divergences associated 

with the hard-scattering kernel can be absorbed into the DA, as discussed in detail in Ref. [46] 



for the case of the pion transition form-factor, leaving us with a result for Iljj n ' which is 
completely finite. Convoluting this renormalised hard-scattering kernel with the asymptotic DA 
results in an expression including terms involving L4 and generalised Nielsen polylogarithms. 
Since we calculate the hard scattering kernel to 0(a 2 Nf), we should take the scale dependence 

of the twist-2 DA to the same order, which involves adding the term 2CfVn(u, v) ln(^ / '/j, )U} ' 

(2) 
to the result for H^ . 

4 Results 

Before coming to our numerical analysis, we must first extract the spectral density from the 
correlation function ILj, and obtain the O(a 2 f3o) QCD sum rules result for the B meson decay 
constant /#. 

4.1 Spectral density 

As in Eq. @ , we define II^ 2 in terms of II^ 2 via 

nj 2 = (PB+p)^l 2 (p 2 B ,q 2 ) + (PB-p)^ T2 (p 2 B ,q 2 ). (35) 

One can then extract the relevant spectral density by taking the imaginary part of the calculated 
correlator, 

p T2 = -Imn^ 2 . (36) 

7T 

An expression for the NNLO correction to px2 is given explicitly in the appendix. As we will 
employ the pole mass for mj in our numerical analysis, we have rewritten the MS mass in terms 
of the pole mass. At 0{a 2 s Nf), this involved adding the term 

A PtI = -CfU^m 2 - 2s) (1(71 + Svr 2 ) + 26 log 4 + 6 log 2 4) ( 37 ) 

O \ £ lit lit J 

to p T ' 2 . Finally, in order to obtain the O(a|/3o) result, Nf in px2 should be replaced by — 3/2/3 - 
Including the contributions at twist-3 to one-loop accuracy and twist-4 to leading order accuracy, 

PU+ (s, 0) = lim (p T2 + /9T3 + Pa + Pp + Pt4 + i ^)' ( 38 ) 

g 2 -s>0 

where /9T3, Pa and p p are contributions at twist-3 and Px 4 are contributions at twist-4 as 
defined in Ref. [16]. An additional twist-4 term, T4 C , cannot be expressed via a dispersion 
relation so must be included separately. Therefore, on taking the Borel transformation of n + , 
we have 

en. 



where we have defined T4c via 



/ dsp u+ (s,0)e- s / M2 +T^\ (39) 

T4( 0) = lim T4 C . (40) 



4.2 Decay constant fs 

Expressing the sum rule as 

f + (p) = -* ( [ S ° dsp n+ (s,0)e( m B- s V M2 +T4 c e</ M2 ] (41) 

we see that a numerical result for /+(0) requires the decay constant fs as input. For consistency 
we use the QCD sum rules result also calculated to O{a 2 s j3o). Although the full 0(a 2 ) corrections 
are sizeable [26l [27] , this is thought to be due to the effect of the classical Coulomb interac- 
tion [28], such that the perturbative expansion is under control. Moreover, the same coulombic 
corrections would also affect the correlator for / + (0)/b- This implies that by employing the sum 
rules result for fs there should be a cancellation between these radiative corrections, as well as 
between the dependence on input parameters such as m& and fi, in / + (0)/# and /#. The QCD 
sum rules result for fs takes the form 

fs= ds p pcr t(s)e^B- s ^ M2 + C qq (qq) + C qGq (qa g Gq)) , (42) 

m B \Jml J 

where C qq and C q c q ar e Wilson coefficients for the operator product expansion (OPE) in terms of 
the quark and mixed condensates respectively |47[ 128] . The spectral density for the perturbative 
contribution p pe it(s) can be expanded in a s , 

P P crt(s) = $U(8) + gp£» + (g) 2 N fP %(8) . . . , (43) 

where the tree level contribution takes the simple form 

Pp e rt ^ = ^-(l-v) 2 - (44) 

The 0(a s ) result /Op ert (s) was obtained from Ref. [48]. The 0{a 2 s ) corrections to p pC rt(s), in the 
case that the light quark is massless, were calculated using Pade approximations and conformal 
mapping and used to obtain semi- numerical results [491 [50] , as an analytical calculation of all 
diagrams was not feasible. We can express /9p er t( s ) m terms of the quantity R F £ s (s), kindly 
provided by the authors of Ref. [59] in publically available code, via 




P f cvt {s) = C F mlsRf^{s). (45) 

To obtain the O(a 2 s fio) result, Nf in /9 per t(s) should be replaced by — 3/2 /3o- The result for 

R F £ s (s) is given at the scale mj, and the pole mass is used for the b quark. We must therefore 
include the O(a 2 f3o) corrections which arise on rescaling a s from m^ to the factorisation scale 
fj,, which take the form 

ApW t (s) = C F ln^$ Tt (s). (46) 

H 1 



Parameter 


Value 


Ref. 


Parameter 


Value 


Ref. 


m n 


139.6 MeV 


[63] 


u 


130.4 MeV 


m 


m B 


5.28 GeV 


[63] 


a s (M z ) 


0.118 


[63] 


m 


0.015 


m 


w 3 


-3 


03] 


m 


10 


m\ 


0J4 


0.2 


03] 


(m) 


(-0.246+0 [J28)3 GeV 3 


m\ 


(qagGq) 


(0.8 ± 0.2) (qq) 


[saES] 



Table 2: Summary of values of parameters used in the numerical analysis. Note the quark 
condensate is given at the scale 1 GeV. 

4.3 Numerical analysis 



From Eq. (J29J) it is clear that making numerical predictions for the twist-2 pion DA comes down 
to determining the Gegenbauer moments. This is only possible via non-perturbative methods 
e.g. QCD sum rules [EEJ EQ [53] or Lattice QCD [H EH [56] . Recently, the UKQCD and RBC 
collaborations computed 02(2 GeV), using Nf = 2 + 1 domain- wall fermions [37]. By combining 
results for a^ (\x) with experimental constraints, i.e. measurements of the 77* it form factor at 
CLEO [35] and CELLO [3J5], an estimate for 04 (//) can be obtained [BD]. However, as this is a 
LCSR calculation, we accordingly adopt 02,4(1 GeV) from Ref. [32] where the LCSR result for the 
pion electro-magnetic form factor [6T] is fitted to experimental data [B2]. The extracted values, 
02(1 GeV) = 0.17±0.08 and 04(1 GeV) = 0.06±0.10, where the errors reflect both experimental 
and theoretical uncertainties, are consistent with other sum rules and Lattice QCD predictions. 
The parameters describing twist-3 and 4 DAs, namely 773, W3, 774 and U4, introduced in Sec. 13.31 
were first calculated in QCD sum rules |52j using non-local operator product expansion and 
conformal expansion. We use the updated results calculated in Ref. [33], as summarised in 
Tab. [21 The error on these parameters is taken to be 50%. The condensates are also required 
as input; we use (qq) and (qcrgGq) as given in Tab. [2] neglecting the gluon condensate as its 
contribution is comparably small. 

Our main numerical analysis is performed using the pole mass nib as input, which we calculate 
to O(a 2 f3o) from the running quark mass. This improves the scale dependence of the final result, 
and avoids any ambiguity in the definition of the lower limit of the integral in Eq. (|4ip . The 
RG improved b quark mass, in the potential subtraction scheme (see Ref. [66]) was calculated 
at NNLO from sum rules in Ref. [67j to be m^ s (2 GeV) = 4.52 ± 0.06 GeV, as in Tab. EJ This 
results in a pole mass of 4.8 GeV at O(a 2 f3o) (and at 0(a 2 )), and in order not to underestimate 
the uncertainty on the pole mass we conservatively adopt mj, = 4.8 ±0.1 GeV. 

The LCSR approach requires a careful choice of numerical values for the continuum limit so 
and the Borel parameter M 2 . We treat the sum rules for /#/ + (0) and Jb separately, obtaining 
independent values of so and M 2 for both. These should be chosen such that the following 
conditions are met: 

• the sum rule exhibits little dependence on, but a clear extremum as a function of these 
parameters; 

• the corresponding sum rule for tub, which can be obtained by differentiating the sum rule 
for f B or /+(0)/ B by 1/Af 2 , is fulfilled to 0.1%, as in Ref. [Ej; 
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0.22 
0.20 
0.18 
0.16 
0.14 



M l 



Figure 2: /b(0) at O(a 2 /3 ) as a function of the Borel parameter M 2 , for central values of input 
parameters (solid) with uncertainties (dotted) calculated as described in the text for /+(0). This 
is compared to the 0(a s ) result calculated using sq = 34.2 GeV 2 (dashed). 



the continuum contribution is under control, i.e. we impose that the integral of the spectral 
density between sq and oo should be approximately 25-30% of the B contribution, between 



in; 



and sq, for / + (0)/b, and 50% for /#; 



as far as possible, the contributions of higher orders in perturbation theory and twists 
should be suppressed. 



Note that we rescale the Borel parameter by (u)" 1 as defined in Ref. [16], as the effective 
Borel parameter in the tree-level sum rule is uM 2 c rather than M 2 C corresponding to M 2 
in Eq. (|41j) . In our numerical analysis we find that so = 34.2 GeV and M 2 = 3.6 GeV for 
Jb, and so = 34.3 GeV 2 and M 2 = 8.1 GeV 2 for /+(0)/b, meet the above requirements. The 
factorisation or renormalisation scale \x is chosen to be the typical virtuality of the b quark, 

as this has previously been found to be an optimal scale |29|. [30l 116] . In Fig. [2] we 



in 



B 



m 



lr 



show fs as a function of M 2 and compare this to the corresponding result at 0{a s ). 
We find that the dominant uncertainties on /+(0) arise due to varying the following: 

• the condensates as indicated in Tab. [21 

• the twist-3 parameter 773 by ±50%; 

• the b quark mass by ±0.1 GeV; 

• the continuum threshold so by ±0.5 GeV and the Borel parameter M2 by ±1.2 GeV for 
both f+(0)f B and f B ; 

• the factorisation scale in the range /x 2 ± 2 GeV . 



The uncertainties arising from each of the above are calculated independently and added in 
quadrature, and we obtain /+(0) = 0.2611 'q 23 . The uncertainties are less than 9%, and could 
be further reduced by better determining the condensates and the twist-3 parameters via, for 
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M' 



Figure 3: /+(0) at O(ajf3o) for central values of input parameters (solid) with uncertainties 
(dotted), compared to the 0(a s ) result calculated using so = 34.3GeV 2 (dashed), as a function 
of the Borel parameter M 2 . 



example, Lattice QCD. Comparing our result for /+(0) to the 0{a s ) result in Fig. [3]shows that, 
despite the ~ 9% 0(a 2 /3o) corrections to fg mentioned earlier, there is little change in /+(0) 
~ 2%. This observation indicates the reliability of the light-cone sum rule approach to the 
calculation of form factors, as it seems that the results are stable with respect to higher order 
corrections. This could further be taken as confirmation that the QCD sum rules result for /# 
should indeed be used in preference to the Lattice QCD result in LCSR calculations of the form 
factors. 

In Ref. |314 [32], /+(0) was calculated using the b quark mass in the MS scheme. Here it 
was argued that this is a natural scheme for the calculation of scattering amplitudes involving 
a virtual b quark at large space- like momentum scales ~ m&. As there are arguments in favour 
of both schemes, we also calculate our result using the MS mass for the b quark. This involves 
replacing the pole mass by the MS mass at the scale fj,, adding NLO corrections found in the 
appendix of Ref. [31] for both twist-2 and 3 scattering kernels. At NNLO translating back to 
the MS scheme for the b mass means removing the correction given in Eq. (137j) . As for fg, 
we take the expressions given in Ref. [26] up to O(a 2 (3o). For the value of the mass, we use 
mb(mb) = 4.19^ ' 06 [68]. Note that, as in Ref. [26], we use the pole mass for the continuum cut- 
off, although using the running mass here instead would change our result negligibly. Imposing 
the same requirements as for the pole-mass scheme, we find so = 35.3 GeV and M 2 = 3.7 GeV 
for f B , and s = 35.7GeV 2 and M 2 = 7.8 GeV 2 for f+(0)f B , and obtain /+(0) 
This is ~ 3% below the result in the pole-mass scheme. 



n oco+0.019 
u - zdz -0.028- 



4.4 Determination of \V u b\ 

As mentioned in the introduction, it is possible to predict \V u b\ using the experimental determi- 
nation of / + (0)|T4f,| and /+(0) from LCSR. In Ref. [21] / + (0)|V^(,| was first obtained by fitting 
various form-factor shape parameterisations to BaBar data [10j . It was observed that the re- 
sults for / + (0)|V^(,| were independent of the parameterisation method chosen. Recently BaBar 
and Belle quote results for /+(0)|V^,&|, extracted by fitting binned data to a Boyd-Grinstein- 
Lebed [17] or Becirevic-Kaidalov [15] parameterisation respectively, as summarised in Tab. [3] 
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Exp. 


No. of bins 


Ref. 


/+(0)|K6| 


\v ub \ 


BaBar 

BaBar 

Belle 


6 
12 
13 


m 

CD 


(1.08±0.06)10- 3 

(8.6±0.3 stat ±0.3 syst )10- 4 

(9.24±0.18 stat ±0.21 syst )10- 4 


(4.13±g;|||th. ±0.23| cxp .)10- 3 
(3-29ig;|| th . ±0.16| cxp .)10- 3 
(3.54+^| th . ±0.11| cxp .)10- 3 



Table 3: Predictions of \V u b\ using / + (0)|V^ft| from analyses in 2010 of B — s> nlv data. 

along with the corresponding value of \V u b\- Where necessary, the statistical and systematic 
uncertainties are added in quadrature. We find that although there is a slight tension be- 
tween Refs m and [12], these predictions are on the whole in keeping with the CKMFitter 



result E3, IK 



>ub\ 



(3.501 



+0.196' 
-0.087. 



10~ 3 and the UTFit result 1701, \V, 



>ub\ 



(3.64 ± 0.11)10' 



They are also in good agreement with the most recent value obtained from LCSR in Ref. [32], 
\V u b\ = (3.50JZ ' 33 |th. =t 0.11|exp.)10 -3 . We look forward to the results from SuperB and Super- 
KEKB which should enable further improvements on the precision of the exclusive determination 
of \V ub \. 



5 Summary 

We have calculated the O(a^/3o) corrections to /+(0) at leading-twist in QCD sum rules on 
the light-cone, and performed a comprehensive numerical analysis of the result, including NLO 
twist-3 and LO twist-4 contributions, leading to a new determination of V u b- We have found 
that in spite of ~ 9% positive NNLO corrections to the QCD sum rules result for fs seen in 
Fig. El the LCSR prediction for /+(0) is stable, increasing by ~ 2% to /+(0) = 0.261+H, 
as shown in Fig. [3j This increases our confidence in the stability of LCSR calculations for 
form factors with respect to this source of radiative corrections, and provides further indication 
that in the calculation of the form factors in LCSR, js should be taken from sum rules rather 
than Lattice QCD. We find that on inclusion of our NNLO correction, the scale dependence 
is reduced, and the main sources of theoretical uncertainty are due to a,2 and rrib- The total 
uncertainty of ~ 9% could be reduced in the future by the determination of the condensates and 
twist-3 parameters on the Lattice. We also evaluate /+(0) using the MS mass for the b quark 
and find /+(0) = 0.252^ ' 02g , in agreement with the result obtained using the pole mass. Finally 
predictions for \V u b\ were obtained in the range (3.29 — 4.13)- 10~ 3 , making use of /+(0)|V^b| from 
BaBar and Belle, in Tab. O We stress that our approach to /+(0) in LCSR is complementary 
to Lattice QCD calculations of f+(q 2 ) as the latter technique is more applicable to the region 
of large q 2 . Therefore the determination of \V u b\ by fitting both our result along with Lattice 
predictions to the combined experimental results [181 [19] would also be of great interest. 
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and the Early Universe" . 

Appendix 

In analogy to Eq. (115|) . we can perturbatively expand the twist- 2 spectral density, 

(0) . a s (i) (u s \ 2 ( 2 ) ,._, 

PT2 = PT2 + ^ PT2 + (^J JV/ A»T2 ■ ■ ■ ■ ( 4 <0 

(2) 

Our NNLO correction p™ then takes the form, 

,3/^,2 



p? ) =/ ,C F {»^llo g 3(l-^) 



/9mj(mj - s) / s \ m b (m^ — s)(20m 4 — 42sm 2 + 7s 2 )\ 2 / to 2 
+ (, s 3 g WJ 27 J log V 1 " V 

/9m!(m 2 -s) 2 / s \ 2to 6 (36to| - 47sm£ + 35s 2 to| - 6s 3 ) /a 

\ s 3 V m2 / ^s 4 V m2 

to^to 2 . - s)(32m| - IOt^sto 2 , - A2Ssm 2 b + 79s 2 ) 



6s 4 



14m 3 (to 2 — s) . /to, 



Li 2 -M log 1 



"°bWi„„(i m b 

s 3 

mliml - s) , , / s \ to?(48to 4 - 134sto? + 107s 2 ) , 

+ ~ Jl Vk lo S — - -^ " THT^ lo S 

3s d \ rn iJ 6s 4 \to^ 

2to 3 (to 2 -s) 1q 2 A_ m| 



s 3 \ s 



_6m fe (TO -s) iQg2 /TO,_ \ +lQg / A 
s d V s / V/^o 

_ / 4to 3 (to 2 - s) /_«_\ + 2to 6 (to 2 -s) \ / _ m|\ 
\ s 3 V m b/ s2 J \ s J 

2to 3 (to 2 - s) 2 / _a_\ _ 2ro| /j_\ + 2to 3 (-15 + tt 2 )(to 2 - s] 



s 3 \ m b/ s2 \ m bJ 3s 3 

2 i „2\ 



2m b {ml - s)(8mf - 8sto^ + s 2 ) / s \ 

s 4 ° g V" 7 ? J 



2m b (6ml + 22sm£ - 19s 2 to^ - s 



+ " v b 1 °- -u 2 

s* \ s 



ml 
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s \ (m b {-2,2ml - 8vr 2 sm 4 + 484sm 4 + 8vr 2 s 2 mg - 239s 2 m 2 . + 24s a 



6 ~0/l OM tfe T lOMI Itj -T O/1 O ,,^-^^0 M6 fe 

mf / v 6s 4 



| 8m 3 (m 2 -s) / s \ | 18m 3 (m 2 -s) L . /m 2 

S 3 \ m 6 / s3 

logf4)f 2m ' (m J" g) log 2 fl-^ 



m£ / V s° \ s 



^/171^ t?\ / O \ \ / TKI^ 



/2mb(m 2 — s)(4m 4 — lOsm 2 + s 2 ) 8m' b J {m b — Sj / ° \ \ i / "% 



\ s 4 s J \ m z 



log ( — ) ) log ( 



2mb(m 2 - s)(7r 2 m 2 - 27m 2 + 9s) 2m 3 , (4m 4 - I4sm 2 + 9s 2 ) 

3s 3 s 4 \ m? 



I2m 3 (™ 2 

^3 



6 v"ift-«J Li / m b 



6 m 3 (m 2 -s) L . / mg\ l0m 3 (m 2 -s) L . /m| 
s 3 V s / s3 V s 

+ ^(l8ir 2 m 6 b - I377r 2 sm£ + 363sm£ - 2l6<(3)m 4 + 98ir 2 s 2 m 2 b 

- 56ls 2 mg + 2l6s 2 C(3)m^ - 3vr 2 s 3 + I98s 3 j I , (48) 

where [1$ is the scale at which the DA moments are calculated. 
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